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ABSTRACT 

Aims. Our aim is to explore the broad-band radio continuum spectrum of LS I +61° 303 during its outbursts by employing the available 
set of secondary focus receivers of the Effelsberg 100 m telescope. 

Methods. The clear periodicity of the system LS I +61°303 allowed observations to be scheduled covering the large radio outburst in 
March-April 2012. We observed LS I +61°303 on 14 consecutive days at 2.6, 4.85, 8.35, 10.45, 14.3, 23, and 32 GHz with a cadence 
of about 12 hours followed by two additional observations several days later. Based on these observations we obtained a total of 24 
quasi-simultaneous broad-band radio spectra. 

Results. During onset, the main flare shows an almost flat broad-band spectrum, most prominently seen on March 27, 2012, where 
- for the first time - a flat spectrum (a = 0.00 ± 0.07, S oc v ff ) is observed up to 32 GHz (9 mm wavelength). The flare decay phase 
shows superimposed “sub-flares” with the spectral index oscillating between -0.4 and -0.1 in a quasi-regular fashion. Finally, the 
spectral index steepens during the decay phase, showing optically thin emission with values a —0.5 to -0.7. 

Conclusions. The radio characteristics of LS I +61°303 compare well with those of the microquasars XTE J1752-223 and Cygnus X- 
3. In these systems the flaring phase is actually also composed of a sequence of outbursts with clearly different spectral characteristics: 
a first outburst with a flat/inverted spectrum followed by a bursting phase of optically thin emission. 
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1. Introduction 

The periodical radio emitting stellar system LS I +61°303 con¬ 
sists of a fast rotating Be star and a compact object of unclear 
nature. One hypothesis is that the compact object is a radio pul¬ 
sar and that the relativistic electrons responsible for the radio 
outburst are produced through an interaction between the rela¬ 
tivistic wind of the pulsar and the wind of the Be star (Maraschi 
|& Treves|1981[|Dhawan et al. |2006[|Dubus|2006| ), i.e. similar to 
the system PSR B1259-63. This system is formed by a 09.5Ve 
star and a compact object with well-detected radio pulses, i.e. 
a radio pulsar. The radio outburst of PSR B1259-63 occur¬ 
ring around periastron passage is explained by electrons accel¬ 


erated at the shock front between the two winds (Chernyakova 
|et al.|[2014 and references therein). For the compact object in 
LS I +61° 303 with the radio outburst towards apastron, how¬ 
ever, an alternative hypothesis suggests an accreting object, i.e. a 
black hole or a neutron star (e.g. Taylor et al. 1992[ |Romero et al. 


|2007| |Massi & Kaufman Ber nado 2009J). In this case, the radio 
emitting electrons belong to a jet typical of microquasars. The 
|Bondi & Hoyle ( |1944] > accretion theory for the eccentric orbit 
of LS I +61°303 with an orbital period P =26.496 days (phase 
<T>) then predicts two accretion maxima and therefore two ejec¬ 
tions of relativistic electrons (|Taylor et al.|1992 Marti & Paredes 


1995 Bosch-Ramon et al. 2006 1 Romero et al. 2001). D uring the 


first ejection near periastron (<E> per j = 0.23 ± 0.02, Casares et al. 

2005} > relativistic electrons suffer severe inverse-Compton (1C) 
losses due to the proximity to the Be star. Consequently, a high 
energy outburst is expected, as indeed confirmed by Fermi/LAV 
observations (JAbdo et al.|2009|), but no or only negligible radio 


emission (|Bosch-Ramon et al.||20061>. For the second accretion 


peak near apastron it has been demonstrated by Bosch-Ramon 
[et al.| ( [2006| ) that the larger distance between the compact object 
and the Be star results in lower IC losses (i.e. a smaller high en¬ 
ergy outburst) and the relativistic electrons can propagate out of 
the orbital plane into a radio jet. Recent Fermi/ LAT observations 
also confirmed this second, smaller high energy outburst (Jaron 
|& Massi|2014| ). 

Synchrotron radiation emitted from one relativistic electron 
population with density N rct gyrating along a magnetic field B 
produces a power-law radio spectrum ( S oc v") with spectral 
slope (Tthin < 0 above a critical frequency Vbreak that depends on 
B and /V rd . Below this critical frequency self-absorption effects 
become important and the emission becomes optically thick with 
a spectral slope cr t hick = 2.5 (e.g. [Kaiser | |2006| . For instance, 
the spectrum of the radio outburst in the pulsar system PSR 
B1259-63 has a spectral index a - -0.7, which is consis¬ 
tent with optically thin synchrotron emission (Fig. 3 in Connors 
et al. 2002|l. In contrast, the jets of microquasars often show 
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flat/inverted radio spectra as demonstrated by Fender ( 2001) >. 
Plasma and magnetic field variations along the jet in fact may 
create regions with different B and N re i values resulting in spec¬ 
tral components with different turn-over frequencies Vbreak- The 
overlap of the optically thin part of one spectral component with 
the self-absorbed part of the adjacent one will result in a flat 
spectrum when observed with a spatial resolution insufficient to 
resolve the jet ( |Kaiser||2006] l. In the microquasar Cygnus X-l, 
radio emission has been measured with a flat spectrum, i.e. a - 
-0.06 + 0.05 and a - 0.07 + 0.04 (Fender et al. 2000| i. In Cygnus 
X-l the black hole is powered by accretion of the stellar wind of 
its supergiant companion star. Since the companion is close to 
filling its Roche lobe, the wind is not symmetric but strongly 
focused towards the black hole (Miskovicova et aT7||201 1 [ and 
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Fig. 1 . Top: Radio outburst of LS I +61°303: flux density (mJy) vs. time (MJD) and orbital phase O as observed with the Effelsberg 
100 m telescope at seven frequencies between 2.6 GHz and 32 GHz. Bottom: Broad-band spectral index a (all frequencies). 


references therein). In contrast, as discussed above in the con¬ 
text of an accretion scenario, accretion in LS I +61° 303 occurs 
in two particular orbital phases due to the eccentric orbit around 
the Be star. Is the accretion phase and subsequent ejection in 
LS I +61°303 able to develop a microquasar jet with a flat spec¬ 
trum? As shown by Corbel et al. ( |2013| >, in the early phase of 
jet formation the low density particles in the jet produce an opti¬ 
cally thin synchrotron power-law spectrum at GHz frequencies. 
Only when the density of the jet plasma increases, a transition to 
higher optical depth occurs resulting in flat/inverted radio spec¬ 
tra (a > 0). In LS I +61°303 strong evidence for a flat radio 
spectrum does exist. Early VLA observations at four epochs re¬ 


vealed a flat spectrum between 1.5 and 22GHz (Gregory et al. 
|1979] >. Furthermore, |Strickman et al.| ( |1998| > carried out multi¬ 
frequency VLA observations sparsely covering one orbit and 
found deviations in the radio spectrum from a simple power law 
during the outburst. Finally, [Massi & Kaufman Bemado ( 2009j l 
measured even inverted, optically thick spectra between 2.2 GHz 
and 8.3 GHz. The primary aim of the current work is to extend 
previous radio observations of LS I +61° 303 and systematically 
study - for the first time - broad-band cm/mm (2.6-32 GHz) ra¬ 
dio spectra and their evolution during a complete outburst. The 
densely sampled observations of LS I +61°303 were performed 
over a period of about three weeks (a total of 24 days) using the 
Effelsberg 100 m telescope at a total of seven frequency bands. 
The paper is structured as follows. In Sect. 2 we briefly present 
the observations and data reduction procedures. In Sect. 3 we 
present the results. Section 4 provides a short discussion and our 
conclusions. 


2. Observations and data reduction 

The multi-frequency flux density measurements of the periodical 
(Jaron & Massij20l3| radio outburst of LS I +61°303 were per¬ 


formed approximately every 12 hours over 14 consecutive days 
between 2012 March 27, orbital phase <t> = 0.31, and 2012 
April 8, <I> = 0.78 (Zimmerman n|2013 i. The orbital phase <I> is 
equal to the fractional part of (t - 43366.275)/7 > or bitai (^orbital = 
26.4960 ± 0.0028 d, Gregory 2002| ). Two additional measure¬ 
ments were obtained on 2012 April 12 (® = 0.94) and 19 (cD = 
0.2) bringing the total observing period to 24 days. At each ob¬ 
serving session lasting for about 2 hours, LS I +61°303 was ob¬ 
served together with primary calibrators consecutively at seven 
frequencies, namely 2.64, 4.85, 8.35, 10.45, 14.60, 23.02, and 
32.00 GHz (110, 60, 36, 28, 20, 13, and 9 mm wavelengths) us¬ 
ing the secondary focus receivers of the Effelsberg 100 m tele¬ 
scope. Flux density measurements were performed using the 
cross-scan technique, i.e. progressively slewing over the source 
position in azimuthal and elevation direction with the number of 
sub-scans matching the source brightness at a given frequency. 
The data reduction from raw telescope data to calibrated flux 
densities/spectra was done in the standard manner as described 
in |Fuhrmann et al.| ( j20081 >. The overall, final measurement uncer¬ 
tainties are of the order of < 1% and <5% at lower and higher 
frequencies, respectively. 


3. Results 

The total flux density evolution at all frequencies is shown in 
Fig. |T] (top). At MJD 56013.2 (0=0.29), the flux density rises 
steeply from around 20 mJy to 100 mJy within 24 hours (un¬ 
til MJD 56014.2) and remains approximately at that level. After 
another 24 hours (MJD 56015.7), it rises yet again to the high¬ 
est flux density level (around 200 mJy at the lowest frequen¬ 
cies) at MJD 56016.7, i.e. the outburst reaches its peak within 
about 3.5 days after onset. We further note that the peak of 
the outburst at frequencies < 10.45 GHz is reached at the same 
epoch, i.e. no significant delay is obvious within our time sam- 
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Fig. 2. Quasi-simultaneous broad-band radio continuum spectra obtained with the Effelsberg 100 m telescope. The spectral index 
a resulting from a power-law fit is indicated in each panel together with the time of observation. The measurement errors are 
comparable to or smaller than the symbol size. We note that up to MJD 56017.2 the spectral index is flat with an average of 
-0.10 ± 0.02 (see text). 


pling of ~ 12 hours. At higher frequencies the situation is less 
clear given the lower data sampling and larger measurement un¬ 
certainties. Subsequently, the slow decay of the outburst takes 
place (MJD > 56017.2) though intersected by several superim¬ 
posed events. These sub-flares occur with rise/decay times of 
~ 12 to 24 hours and with amplitudes of a few lOmJy. We note 
that during the last sub-flare (around MJD 56024), the peak at 
higher frequencies (here: 10.45 and 14.60 GHz) appears to oc¬ 
cur 24 hours earlier. In Fig. [2] the broad-band radio spectra of 
all epochs are shown. In order to quantify the spectral behaviour 
of LS I +61°303, broad-band spectral indices a (S oc v") were 
obtained by performing a power-law fit to all available multi¬ 
frequency flux densities of a given epoch and the resulting fits 
are superimposed in Fig. [2] The spectral index evolution with 
time is shown in Fig. [T| (bottom: a vs. MJD) and is character¬ 
ized by an overall long-term decreasing trend with (i) the flattest 
spectral indices occurring during the early raise/peak of the out¬ 
burst and (ii) increasingly steeper values during the outburst de¬ 
cay with the steepest values reached towards the end (a —0.5 


to -0.7). A linear fit to the spectral indices reveals that on av¬ 
erage a spectral steepening of Aa^O.03 per day takes place 
over the total observing period of 24 days. Examining the spec¬ 
tral index behaviour in more detail we basically identify three 
stages (vertical dotted lines in Fig. |T] bottom). During the rise 
until shortly after the main peak, i.e. from MJD 56013.2 to MJD 
56017.2, the broad-band spectrum remains almost flat (average 
spectral index a - -0.10+0.02) and an even inverted spec¬ 
trum is observed at high frequencies during some epochs (i.e. 
0+4 GHz/23 GHz ~0.6, 0.4, and 0.6 at MJD 56013.8, MJD 56015.2, 
and MJD56015.7, respectively). Subsequently, a stage covering 
the decay phase follows with the superimposed sub-flares lasting 
until MJD 56023.1 and with spectral indices dropping to values 
as low as -0.39 (mean value: -0.29+0.03). Here, however, a ap¬ 
pears to oscillate between -0.4 and -0.1 in a quasi-regular fash¬ 
ion, whereas the local peaks in spectral index roughly coincide 
with the peaks of the sub-flares seen in the total intensity light 
curves (Fig. [T] top). We further note that also during the decay 
stage several epochs show an inverted spectrum at high frequen- 
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cies. During the last stage of the decay (MJD 56025 to MJD 
56037), the spectral index finally drops to the steepest observed 
values ranging between -0.5 and -0.7, i.e. clearly optically thin, 
with an average of -0.5 ±0.1 during this period. 


4. Discussion and conclusions 

During our Effelsberg campaign in 2012 March-April, 
LS I +61°303 was monitored at seven frequencies between 2.6 
and 32 GHz over a total period of 24 days with a cadence of 
about 12 hours (over the first 12 days). Our observations reveal 
the following results: 


1 . 


The broad-band data show a main peak with nearly flat spec¬ 
tra followed by a sequence of small amplitude, optically thin 
events. This sequence of outbursts with different spectral 
characteristics is also found in Cygnus X-3: Fig. [3|r shows 


GBI observations of Cygnus X-3 (Waltman et al. 1995 1 , 


2 . 


where a first outburst with almost flat spectrum (a ~ 0.05) 
is seen, followed by a minor optically thin flare. Figure 
[3J: displays a first outburst exhibiting an inverted spectrum 
(a ~ 0.6) followed by a bursting phase of optically thin emis¬ 
sion, as also seen for LS I +61°303 in the GBI observations 
of Fig. [3}t. Clearly both sources, LS I +61°303 and Cygnus 
X-3, show the same sequence of events: a first outburst with 
flat/inverted spectrum is followed by one or more optically 
thin flares. A third microquasar with the same complex out¬ 
burst behaviour is XTE J1752-223 ( jBrocksopp et al.[2013) l. 
In Fig. [3^1 we show one outburst at 5.5/9 GHz (see also Fig. 1 
in Brocksopp et al .|2013[ >, where the first peak with a ~ 0.2 
is clearly distinguished from the subsequent optically thin 
flares. Brocksopp et al. (2013 1 associate the first flare with 
the compact XTE J1752-223 jet and the sequence of opti¬ 
cally thin flares with discrete ejections or shocks. 

During the sequence of LS I +61°303 sub-flares presented in 
Fig. E a appears to oscillate in a quasi-regular fashion. Such 
spectral variability has so far only been discussed in the mi¬ 
croquasar GRS 1915+105 (see Fig. 2 in Fender et al. |2002| l. 
As apparent from the inset of Fig. [3je (2.2/8.3 GHz spectral 
index vs. time), Cygnus X-3 also shows similar spectral in¬ 
dex oscillations during the three, minor sub-flares occurring 
between MJD 48470 and MJD 48480. 


The dashed line in Fig. [3jt indicates the epoch where [Paredes | 


et al. (2000) obtained a steep LS I +61°303 cm/mm spec¬ 


trum up to 250 GHz (1998 March 15) during one of the opti¬ 
cally thin sub-flares following the main flare. As a follow-up 
of our work, it is certainly of interest to observe the main out¬ 
burst of LS I +61°303 also at higher radio frequencies. Indeed, 
flat/inverted spectra in microquasars have also been established 
at short millimeter wavelengths (jFender et al. 2000| Fender 


2001) and it will be important to clarify if the short accretion 


phase in the eccentric orbit of LS I +61° 303 is able to maintain 
a flat spectrum up to such high frequencies. 
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